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SUMMARY 


A wind tunnel experiment on the 10-percent-scale model of the F-16XL aircraft included 
longitudinal static tests, and oscillatory and ramp tests in pitch. Static tests investigated the effect 
of angle of attack, slideslip angle and control surface deflection on aerodynamic coefficients. For 
dynamic testing and analysis, the report presents data from small-amplitude oscillatory tests at 
nominal values of angle of attack between 20° and 60°, five frequencies from 0.6 to 2.9 Hz and 
one amplitude of 5°. A simple harmonic analysis provided Fourier coefficients associated with the 
in-phase and out-of-phase components of the aerodynamic coefficients. A strong dependence of 
wind tunnel oscillatory data on frequency led to the development of models with unsteady 
aerodynamic terms in the form of indicial functions. Two models expressing the variation of the 
in-phase and out-of-phase components with angle of attack and frequency were proposed and their 
parameters estimated from measured data. Both models with the estimated parameters were in 
good agreement and represented unsteady effects observed in the aerodynamic forces quite well. 
The estimated parameters were close to the results from the static test and both models 
demonstrated good prediction capability. 
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SYMBOLS 


A 

a,b t ,c 

a,b 
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C La ( t\C Lq (r) 
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f 
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k 

£ 

m 

n 

n c 

q 

s 

S 
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Tx 
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u,v 


amplitude ratio 

parameters in indicial function 
Fourier coefficients (Appendix A) 
wing span, m 

lift, normal force, pitching-moment coefficients 
indicia] functions 

wing mean aerodynamic chord, m 
frequency, Hz 
cost function 

reduced frequency, k = c o£ IV 
characteristic length, i = c / 2 
number of frequencies 
number of angles of attack 
number of cycles 
pitching velocity, rad/sec 
standard error 
wing area, m 2 
period, sec 

time constant, 7] =X\(/V, sec 
time, sec 

input (Appendix A) 

aerodynamic derivatives in model for in-phase and out-of-phase 
component of oscillatory data, respectively 
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u,v 

V 
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Z u »^v 

a 

0 

A 

Se 

0 

X 

'Cl 

CO 


in-phase and out-of-phase components of oscillatory data, 
respectively 

airspeed, m/sec 

output (Appendix A) 

terms in model equation (10) and (11), respectively 
angle of attack, rad or deg 
sideslip angle, rad or deg 
increment 

elevon/flaperon deflection, rad or deg 

parameter in Model II 
time delay, sec 

nondimensional time constant, Xj = V/ tb\ 
phase angle, rad or deg 
angular frequency, rad/sec 


Superscript: 


estimated value 


Superscript over aerodynamic derivative: 

oscillatory data 

Subscript: 

A amplitude 

0 nominal value 


Aerodynamic derivatives: 




for A = L,N, or m 


qc qc 2 ac „ 
a = — , — T ,a , — , or o„ 
2V 4V 2 2V e 
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INTRODUCTION 


For better understanding of aircraft aerodynamics in large amplitude maneuvers, NASA 
Langley Research Center conducted a series of wind tunnel tests on a model of the F-16XL 
aircraft. These tests included measurement of aerodynamic forces and moments under static 
conditions followed by oscillatory and ramp tests in pitch. Static measurements allowed 
investigation of changes in aerodynamic coefficients due to angle of attack, sideslip and control 
surface deflection. Dynamic testing quantified changes in aerodynamic coefficients with frequency 
and amplitude. 

The purpose of this report is to 

(a) document static and small-amplitude oscillatory data for qualitative assessment; 

(b) determine a mathematical model with unsteady aerodynamic terms for oscillatory data at 
different angles of attack and frequencies. 

The parameters in a model for small-amplitude oscillatory data were estimated by a nonlinear 
estimation technique of reference 1. Then, a new procedure using splines for model formulation 
and linear regression was applied. 

The report begins with a description of the experiment. The results from static and oscillatory 
test are then described followed by model structure determination and parameter estimation from 
oscillatory data at various angles of attack and frequencies. The report is completed by concluding 
remarks. 

MODEL AND TESTS 

A three-view sketch of the 0.10-scale F-16XL model is shown in figure 1 together with some 
basic model dimensions. Static and dynamic tests were conducted in the NASA Langley 
12-Foot Low-Speed Wind Tunnel. For both tests, the model was mounted on a dynamic test rig 
through a six-component strain-gauge balance. The dynamic test rig is a computer controlled, 
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hydraulically-actuated system which was sting-mounted on a C-strut support system. The 
mounting arrangement rotated the model about the reference center of gravity location of 0.558 c , 
over an angle of attack range of -5° to 80°. The maximum capability of the dynamic test rig was 
260 deg/sec pitch rate and 2290 deg/sec 2 pitch acceleration. Further description of the dynamic test 
rig may be found in reference 2. The tests were conducted at a dynamic pressure of 192 Pa (4 psf) 
resulting in a Reynolds number of 10 6 based on the mean aerodynamic chord. 

All data were obtained with the leading-edge flaps at 0° deflection. Static data were obtained 
for angles of attack from -5° to 80° at zero sideslip and zero deflection of trading-edges surfaces 
(flaperons and elevons). Additional data were obtained for four values of sideslip, ((3 = -5°, 
—10°,— 20°, and -30°), and for two values of symmetric trading-edge control surface deflections 
(8 e = -20° and 20°). Oscillatory data were obtained at initial angles of attack (do) between 20° 
and 60° at five frequencies (0.6, 1.0, 1.41, 1.75, and 2.94 Hz) and an amplitude of ±5°. The 
effects of sideslip and control surface deflections were measured only at an initial angle of attack of 
35° and two frequencies (0.6 and 1.41 Hz) at the same sideslip values and control deflections as in 
the static test. 

Data recorded during the test runs included a linear variable differential transformer reading 
from the dynamic test rig to determine the pitch angle, six-component force and moment data from 
a strain-gauge balance, and wind tunnel dynamic pressure. Data were sampled at 100 Hz with an 
in-line 100 Hz low pass filter. All data channels were subsequently filtered using a 6 Hz low pass 
filter. 


STATIC DATA 


Among the measured and computed data only three aerodynamic coefficients, C L ,C N , and 
C m , are included. These coefficients are plotted against the angle of attack in figure 2. The 
variation of C L with a is almost linear for a < 25°, the maximum lift occurs at a = 34°. The 
normal force is also changing linearly with a till a = 25°, and reaches maximum value at 
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a ~ 38°. The model exhibits static instability for a < 30°. This pitch-up tendency is inherent to 
highly swept wings at low speeds as documented in many publications, e.g., reference 3. 

Figure 3 shows the effect of sideslip on longitudinal characteristics. For J3 < 5° there arc 
almost no changes in these characteristics. Larger sideslips, however, cause a decrease in values 
of all coefficients, mostly in the region of a between 20° to 40°. Pitching moment decreases due 
to sideslip at all angles of attack > 20°. Large sideslip angles eliminate the pitch-up tendency in 
small segments of the angle of attack in the pre-stall region. The effect of symmetric elevon and 
flaperon deflection is presented in figure 4. In the post-stall region the control effectiveness is 
larger for negative (trailing edge up) deflection than for positive. The stability and control 
derivatives evaluated from the data in figure 2 and 4 with zero sideslip angle are plotted in figure 5 
and 6. The control derivatives were computed for both positive and negative control deflections. 
From these plots a rapid decrease in control effectiveness with increased a in the post-stall region 
is apparent. All the measured and computed data in figures 2 to 6 are summarized in tables I to IV. 

OSCILLATORY DATA 

As indicated earlier, the basic set of oscillatory data was obtained at five frequencies, nine mean 
values of the angle of attack, a 0 , and one amplitude, = 5°. The following table show the 
values of frequencies and the corresponding reduced frequencies selected for the test. 


to, 

rad/sec 

f, 

Hz 

, t 
k = —to 

V 

3.77 

.60 

.081 

6.28 

1.00 

.134 

8.86 

1.41 

.190 

11.00 

1.75 

.237 

18.47 

2.94 

.397 
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The mean values of a were 20.8, 25.9, 30.8, 35.8, 40.8, 45.8, 50.8 55.9 and 61.1 with the 
maximum deviation from these values during the repeated settings less than 0 . 1 °. 

The examples of the oscillatory data at one frequency and two mean angles of attack are given 
in figure 7 and 8 . For a 0 = 30.8°, the plot of C L (t) indicates a small change from a sine wave, 
whereas the plot of C m (t) shows the existence of at least one additional harmonic. For a 0 = 50.8° 
the departure of C L (t) from a sine wave is more pronounced and C m (t) does not resemble a 
simple harmonic motion at all. Further examination of the remaining time histories revealed severe 
distortion of measured pitching-moment coefficient at a > 40° at all frequencies, and the response 
can be characterized as a chaotic. The reason for this behavior has not been found yet. 

For the following analysis of the oscillatory data it is assumed that the aerodynamic coefficients 
are linear functions of the angle of attack, pitching velocity and their rates. Then, for example, the 
increment in the lift coefficient with respect to its mean value can be formulated as 



where for the harmonic motion 

Aoc = sin cot 
a = <7 = toot ,4 cos cor 
d = q = -Gi 2 a A sin cor 

Substitution of (2) into (1) yields 

ACi = — k 2 Ci jsincor + a^&lo,. + C^ jcoscor 

= o. A {c La sincor + £C^ cos cor j 


( 2 ) 


( 3 ) 
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where 


a A C L a ~ a A( C L a _ * 2 c £vj) 
a A kC L q =a A k {^L 6i +C L q ) 


(4) 

(5) 


represent the Fourier coefficients. The in-phase and out-of-phase components of and 

Q, ) can be obtained by integrating the time histories of AC L over n c cycles as 


Cr =• 


2 "c T 




|AC^(Osinw/dt 


C, = 


2 


( 6 ) 


^ a A k/i c r 0 


jAC L (r)coscot d/ 


where T = 2k/co . In the present analysis the integrals in equation ( 6 ) were obtained by a 
numerical technique of reference 4 applied to the last three cycles of the steady harmonic outputs. 

The in-phase and out-of-phase components of all three coefficients are presented in tables V to 
VII and plotted against the angle of attack in figure 9. The figure shows the effect of frequency 
which is especially strong on the out-of-phase components. The variation of these components 
with the reduced frequency for two values of (Xq can be seen in figure 10. The plot of C m ^ for 
ct 0 - 50.8° may indicate the problem with measured C m (t) mentioned earlier. 

Figure 1 1 and 12 compare the measured coefficients C L and C m at one frequency and two 
values of (Xo with those computed from eq. (3) and (4). These coefficients are plotted against the 
angle of attack with time as a parameter. As shown in Appendix A, these plots should be ellipses 
if the formulated model given by eq. (1) is correct. Any deviation of measured C L (a,t) and 
C m (a,t) from this shape might be caused by aerodynamic nonlinearities and/or flow separation. 
Substantial differences between measured and computed pitching-moment coefficient exist as 
indicated by results in figure 1 1 and 12. When a and C L are considered as the input and output 
variables of a dynamical system then their relationship is steady harmonic motion can be defined by 
the amplitude ratio and phase angle, see Appendix A. The frequency response of these quantities 
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for two values of otp is presented in figure 13. From these plots a strong effect of frequeny on the 
amplitude and phase within the range of k from 0 to 0.2 is apparent. 

A comparison of steady and oscillatory data for all three coefficients is given in figure 14. The 
oscillatory data presented are for one frequency and five values of ar 0 from 20 to 60°. For 
a 0 = 20° the reduction of ellipses C L ( a) and C N ( a) to almost a line is caused by small values of 
and C N q . Possible reasons for distorted ellipses has already been discussed. The effect of 
sideslip, elevon/flaperon deflection and tail-off configuration on the components of all three 
coefficients are summarized in figure 15 and 16. 

PARAMETERS ESTIMATED FROM OSCILLATORY DATA 

A strong dependence of wind tunnel oscillatory data on frequency lead to the development of 
models with unsteady aerodynamic terms, see e.g. reference 5 and 6. As an example, the model 
for the lift increment was formulated in reference 4 as 

= } C La (/ - x) j- adx + £ } C La q(x)dx (7) 

where (f) and Ci ( t ) are the indicial functions. 

For obtaining the model with limited number of parameters, it was assumed that the effect of 
q(t) on the lift could be neglected and the indicial function C L ^ ( t ) can be expressed as 

C La (t) = a(l-e- b '‘) + c (8) 

Considering the above mentioned assumptions, equation (7) is simplified as 

AC l =C l ^-a(x)dx + ^-C L {~>)q{t) (9) 

“ dx V q 

where (®) and Ci (°°) are the rates of change of Q, with a and q in steady flow. The 
steady form of equation (9) for harmonic changes in a is identical to that of equation (3), that is 
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AC L (t) = a A^L a sin©/ + o. A kC L ^ cos©/ 


( 3 ) 


However, with the indicial function of equation (8), the expressions for C L ^ and C L ^ have the 
form 



= C L a H-fl 


x?* 2 

1 + X 2 ifc 2 


( 10 ) 



= C L q ( 00 )- a 


l + xf/t 2 


( 11 ) 


where Tj = Vjb\i is the nondimensional time constant. 

From the experiment, the in-phase and out-of-phase components are usually obtained for 
different values of the angle of attack and reduced frequency while keeping the amplitude of the 
oscillations constant. Then the equations (10) and (11) can be generalized as 


Uji =«,• -a iZu . 


Vji=v i -a i z v . 


( 12 ) 

(13) 


These equations define Model I where for the lift coefficient 


u i=CL a < a i) 

T?*? 


v i= c L t < a i) 


1 + X 


2 1 2 
1 k j 


Zv, = 


It 


i+X^ kj 


for 


i = l,2,...,n 
j = l,2,...,m 


In equations (12) and (13) there are, in general, 3n + l unknown parameters: and x j . 

They can be estimated from experimental data Uji and Vji by minimizing the cost function 

m n [- 

Ji = x x 

y=i i=i l L 
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More about the estimation procedure can be found in reference 1 . 

In formulating airplane equations of motion it might be more convenient to obtain expressions 
for u, v, and a as a function of the angle of attack rather than their discrete values. For that reason 
the previous model was reformulated as Model II defined by the following equations: 


JF 

i 

JO, 

'IS' 

II 

(15) 

N? 

jjf 

"o' 

i 

jiT 

ii 

■ 

(16) 


The form of expressions for «(a,), v(a, )and a(a t ) can be either specified from the variation of 
estimated parameters in Model I or the form can be postulated in terms of polynomials and/or 
polynomial splines. In the second case adequate models for u(a),v(oc) and a(a) can be 
determined from measured data by a stepwise regression. The previously estimated value of 
parameter Tj can be used as an a priori value thus simplifying the parameter estimation procedure 
using Model IT and the cost function 


m fl fr .-12 

J II = X Xj u ji — ^u(otj) — a(ot,)z,y)J 


j=U=\ 


h [ v ji -(v(a { )-fl(a,)z w )] 2 } 


(17) 


Both estimation procedures were applied to oscillatory data at four frequencies and nine value 
of the angle of attack. The measured data for frequency k = .190 were not included in the 
estimation. This set of data was used to demonstrate the model prediction capabilities. The 
estimated parameters and in Model I are plotted in figures 17 to 19. In these figures, the 

2s confidence intervals for parameter estimates are also included. The minimum and maximum 
values of standard errors are summarized in Table VTTI. The parameters C La (°°) and C Na (°°) are 
in good agreement with those from static data displayed in figure 5. There is, however, large 
discrepancy between C m ^ (°°) from oscillatory and static data caused by problems in measured 
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time histories C m (t), as discussed earlier. No values, theoretical or experimental, were available 
for derivatives Cl q ’ C N q and C m q for comparison with results from oscillatory tests. The 
parameter values a,- indicate smooth variation of unsteady terms with a and the largest effect of 
unsteady terms on the coefficients at a around 40°. The estimated values of Tj , and computed 
values of parameter b\ and 7] are presented in Table IX together with their standard errors. These 
results indicate that the time constant associated with the unsteady effects is about 0.4 sec. for the 
lift and normal force, and about 0.6 sec for the pitching moment. 

A comparison of measured in-phase and out-of-phase components of all three coefficients with 
those estimated is given in figures 20 to 22. The predicted components computed from model 
equations (10) and (11) for £ = .190 are presented in figure 23 to 25 together with the 
corresponding measured values. These figures demonstrate that Model I for C L and C N is a good 
predictor. No conclusion about predicted values of C m ^ and C m ^ can be made because of chaotic 
behavior of measured values C m (t) at a > 40°. 

Expressions for terms u(a),v(a) and a(a) in Model II for the lift coefficient were selected 
from results in figure 26 as 

C L(x (a) = «(a) = 0 o +0]a + 02a 2 +03(a-O.8O3) 2 

C l (a) = v(a) = 0 4 +0 5 a + 0 6 a 2 (18) 

a(a) = 0 7 +0 8 a + 09a 2 + 0 lo (a-O.8O3) 2 

where 

(a - 0.803) 2 = 0 for a < 0.803 

= (a-0.803) 2 for a >0.803 

An identical model was postulated for the normal-force coefficient. There was no attempt to 
estimate parameters in the pitching-moment equations. 

The parameters in equations (18) were obtained by a linear regression for given values of 
parameter X| . The mean values and standard errors of the estimates are summarized in table X. 
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The results indicate parameter uncertainty in expressions for (a) and C Ng (a). Low accuracy 
of these parameters could result from small number of data points and their small sensitivity due to 
small magnitudes of and C N q • The estimated terms of the in-phase and out-of-phase 
components are plotted in figures 26 and 27, and compare with estimated parameters in Model I. 
The agreement between both sets of results is very good. The measured and estimated 
components of C L and C N are presented in figure 28 and 29. With the exception of Cl u and 
C N a for k = 0.397 and 0.081, the estimated components are close to the measured values. 
Finally, figure 30 and 31 demonstrate prediction capabilities of Model II which are as good as 
those of Model I. 


CONCLUDING REMARKS 

A wind tunnel experiment on a 10-percent-scale model of the F-16XL aircraft was conducted at 
NASA Langley Research Center. This experiment concluded longitudinal static and oscillatory 
tests. Static tests investigated the effect of the angle of attack, sideslip angle and control surface 
deflection on aerodynamic coefficients from which only the lift, normal-force and pitching moment 
are presented in this report. The variation of lift with angle of attack was found to be almost linear 
for angles of attack less than 25°; the maximum lift occurred at an angle of attack of about 34°. 
Similar results were obtained for the normal force. The model exhibited static instability for angles 
of attack less than 30°. For small sideslip angles, less than 5°, there were almost no changes in the 
longitudinal characteristics. Large sideslips, however, caused a decrease in values of all 
coefficients. Changes in symmetric elevon and flaperon deflections indicated that in the post-stall 
region the control effectiveness was larger for negative (trailing-edge-up) deflections than for the 
opposite defections. 

For dynamic testing and analysis, the report presents only data from small-amplitude 
oscillatory tests at nominal values of angle of attack between 20° and 60°, five frequencies from 
0.6 to 20.9 Hz and one amplitude of 5°. A limited amount of data exhibiting the effect of sideslip 
angle, control surface defection and tail-off configuration is also included without any discussion. 
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Time histories of the lift and normal force coefficients for single of attack greater than 30° indicated 
relatively large changes from a sine wave, probably caused by aerodynamic nonlinearities. The 
time histories of the pitching moment, however, did not resemble a simple haromonic response at 
all. With increasing angles of attack these time histories were getting closer to chaotic motion. 
Regardless of these irregularities in time histories, it was assumed that for data analysis the 
aerodynamic coefficients were linear functions of the angle of attack, pitching velocity and their 
rates. Then, a simple harmonic analysis over three cycles of oscillation provided the Fourier 
coefficients associated with the in-phase and out-of-phase components of the coefficients. These 
components were plotted against the angle of attack for a given value of reduced frequency. 

A strong dependence of wind tunnel oscillatory data on frequency led to the development of 
models with unsteady aerodynamic terms in the form of indicial functions. These functions were 
postulated as simple exponentials where the unknown parameters included aerodynamic 
derivatives, the exponents and multiplication terms. Two models expressing the variation of the 
in-phase and out-of-phase components with the angle of attack and frequency were proposed. In 
the first model, the unknown parameters, including the exponent, were estimated from measured 
data by a nonlinear least squares technique at nominal values of angle of attack. In the second 
model, polynomial splines expressing the variation of unknown parameters with the angle of attack 
were formulated and their parameters estimated by linear least squares technique for a known value 
of the exponent. Both models with the estimated parameters were in good agreement and 
represented unsteady effects observed in the lift and normal force quite well. In addition, some of 
the estimated parameters were close to the results from the static test and both models demonstrated 
good prediction capabilities. 
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APPENDIX A 


HARMONIC MOTION 

Let the input to a linear, time-invariant, dynamical system be 

u = a 4 sin (at 

Then the steady state output of the system has the form 

y = asin (at + bcos(at 
= Asin(tot-(j)) 

where 

a = Acos<() 
b = -Asin<j) 

and 

A 2 = a 2 +b 2 
b 

tan<(> = — 
a 

Replacing y by AC^ and u by a it follows from equations (A3), (A4) and (3) that 

“ = u A C La 

b - a A kC L g 

and the expression for amplitude ratio and phase angle takes the form 

A = (X A 

<J> = tan' 



(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 
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Returning to equations (Al) and (A2), and considering unit amplitude of the input variable 

y = a sin cot + b cos to/ 

= a sin tot + b4 1 - sin 2 tot 
= au + by 1 - w 2 

Taking square of both sides of (A8) yields 


(A8) 


y 2 - 2auy + ^ a 2 + b 2 jw 2 = b 2 


(A9) 


Replacing again y by A C L and u by a in equation (A9), and using, equation (A5) gives 


A C 2 h - 2C l n AC l cl + + k 2 C 2 ^ jtx 2 = 


k 2 C 2 


L q 


(A10) 


The discriminant of eq. (A 10) is 


-C 


~Cr 


L a 

ct a +k 2 c 2 


= k 2 Cl^ >0 


— 2 

which means that eq. (A10) represents an ellipse. For kC ^ = 0 the ellipse is reduced to a line as 
follows from (A8). 
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0.7470 

0.7447 

0.7181 

0.6625 

0.61 19 

1.0309 

0.3637 

20.0000 

0.9774 

0.9831 

0.9672 

0.8482 

0.8041 

1.2382 

0.6066 

22.0000 

1.0543 

1.0919 

1.0545 

0.8800 

0.8358 

1.3285 

0.6879 

24.0000 

1.1468 

1 .1593 

1.1202 

0.9199 

0.8523 

1.3994 

0.7718 

26.0000 

1.2070 

1 .2258 

1.1604 

0.9724 

0.8629 

1.4422 

0.8462 

28.0000 

1.2529 

1.2716 

1.1700 

1.0138 

0.8374 

1.4577 

0.9033 

30.0000 

1.2726 

1 .2718 

1.1916 

1.0029 

0.8447 

1.4314 

0.9484 

32.0000 

1.2780 

1.2696 

1.1907 

0.9668 

0.8594 

1.3981 

0.9811 

34.0000 

1.2650 

1.2613 

1.1988 

0.9926 

0.8588 

1.3553 

1.0126 

36.0000 

1.2388 

1.2233 

1.1806 

1.0152 

0.8556 

1.3389 

1.0192 

38.0000 

1.2143 

1.1374 

1.1376 

1.0244 

0.8682 

1.2628 

1.0033 

40.0000 

1.1347 

1.0445 

1.0513 

1.0237 

0.8860 

1.1 145 

0.9366 

45.0000 

0.9255 

0.8813 

0.8833 

0.8518 

0.8402 

0.9767 

0.7748 

50.0000 

0.8377 

0.8235 

0.8153 

0.7826 

0.7547 

0.8798 

0.7037 

60.0000 

0.7202 

0.7046 

0.6852 

0.6658 

0.6289 

0.7027 

0.5963 

70.0000 

0.5403 

0.5240 

0.5048 

0.4967 

0.4616 

0.5004 

0.4553 

80.0000 

0.2950 

0.2894 

0.2850 

0.2742 

0.2567 

0.2331 

0.2629 


Table I. Effect of angle of attack, sideslip, and elevons and flaperons on lift coefficient. 
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a, deg 


P = - 30 - 
5 - = 0 - 


- 4.0000 

- 0.1022 

- 0.1075 

- 0.1064 

- 0.1123 

- 0.0899 

0.1980 

- 0.4403 

0.0000 

0.0535 

0.0500 

0.0478 

0.0380 

0.0427 

0.3555 

- 0.2586 

5.0000 

0.2612 

0.2568 

0.2505 

0.2358 

0.2274 

0.6055 

- 0.0639 

10.0000 

0.5058 

0.5040 

0.4858 

0.4499 

0.4283 

0.8426 

0.1469 

15.0000 

0.7714 

0.7690 

0.7414 

0.6847 

0.6346 

1.0769 

0.3875 

20.0000 

1.0335 

1.0400 

1.0232 

0.8991 

0.8545 

1.3297 

0.6553 

22.0000 

1.1269 

1.1689 

1.1291 

0.9456 

0.9008 

1.4459 

0.7512 

24.0000 

1.2416 

1.2575 

1.2155 

1.0035 

0.9333 

1.5457 

0.8534 

26.0000 

1.3256 

1.3495 

1.2794 

1.0775 

0.9608 

1.6195 

0.9491 

28.0000 

1.3993 

1.4245 

1.3135 

1.1430 

0.9508 

1.6677 

1.0306 

30.0000 

1.4480 

1 .4526 

1.3638 

1.1467 

0.9779 

1.6715 

1.1021 

32.0000 

1.4839 

1.4809 

1.3915 

1.1362 

1.0147 

1.6692 

1.1637 

34.0000 

1.5012 

1.5044 

1.4320 

1.1921 

1.0360 

1.6567 

1.2267 

36.0000 

1.5046 

1.4968 

1.4447 

1.2478 

1.0557 

1.6773 

1.2634 

38.0000 

1.5115 

1.4311 

1.4298 

1.2908 

1.0972 

1.6268 

1.2749 

40.0000 

1.4546 

1.3533 

1.3597 

1.3254 

1.1484 

1.4834 

1.2268 

45.0000 

1.2840 

1.2354 

1.2374 

1.1949 

1.1730 

1.4104 

1.1051 

50.0000 

1.2674 

1.2604 

1.2460 

1.2033 

1.1561 

1.3962 

1.0979 

60.0000 

1.3642 

1.3566 

1.3196 

1.2935 

1.2230 

1.4325 

1.1702 

70.0000 

1.4411 

1.4347 

1.3916 

1.3736 

1.2955 

1.5044 

1.2660 

80.0000 

1.4643 

1.5223 

1.4958 

1.4383 

1.3444 

1.5504 

1.3743 


Table II. Effect of angle of attack, sideslip, and elevons and flaperons on normal-force 

coefficient. 
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a , deg 

p = o- 

fie =0- 

P = - 5 * 
8 e = 0 ‘ 

P = - 10 ' 

fie = 0 " 

P = - 20 ' 
fi 8 = 0 - 

P = - 30 - 

fie = 0 ' 

P = O ' 
8 e = 20 ' 


- 4.0000 

0.0201 

0.0156 

0.0265 

0.0163 

0.0237 

- 0.0627 

0.0937 

0.0000 

0.0242 

0.0270 

0.0318 

0.0261 

0.0329 

- 0.0558 

0.0998 

5.0000 

0.0441 

0.0458 

0.0476 

0.0428 

0.0515 

- 0.0479 

0.1204 

10.0000 

0.0733 

0.0744 

0.0755 

0.0702 

0.0764 

- 0.0095 

0.1525 

15.0000 

0.1107 

0.1107 

0.1 110 

0.1040 

0.1035 

0.0385 

0.1940 

20.0000 

0.1530 

0.1531 

0.1515 

0.1458 

0.1051 

0.0855 

0.2359 

22.0000 

0.1685 

0.1734 

0.1679 

0.1605 

0.1248 

0.1047 

0.2521 

24.0000 

0.1858 

0.1880 

0.1841 

0.1574 

0.1470 

0.1246 

0.2694 

26.0000 

0.2009 

0.2059 

0.1970 

0.1726 

0.1627 

0.1437 

0.2834 

28.0000 

0.2161 

0.2214 

0.1995 

0.1938 

0.1577 

0.1603 

0.2931 

30.0000 

0.2244 

0.2243 

0.2094 

0.1985 

0.1476 

0.1680 

0.2977 

32.0000 

0.2256 

0.2238 

0.2168 

0.1771 

0.1377 

0.1700 

0.2966 

34.0000 

0.2214 

0.2222 

0.2196 

0.1738 

0.1357 

0.1708 

0.2987 

36.0000 

0.2183 

0.2164 

0.2162 

0.1822 

0.1407 

0.1773 

0.2955 

38.0000 

0.2186 

0.2020 

0.2068 

0.1863 

0.1460 

0.1698 

0.2891 

40.0000 

0.1995 

0.1890 

0.1943 

0.1926 

0.1479 

0.1478 

0.2631 

45.0000 

0.1867 

0.1728 

0.1750 

0.1701 

0.1664 

0.1594 

0.2431 

50.0000 

0.1793 

0.1634 

0.1596 

0.1462 

0.1386 

0.1596 

0.2386 

60.0000 

0.1578 

0.1493 

0.1386 

0.1371 

0.1144 

0.1374 

0.21 14 

70.0000 

0.1208 

0.1181 

0.1065 

0.1074 

0.0794 

0.1093 

0.1706 

80.0000 

0.0594 

0.0466 

0.0428 

0.0560 

0.0479 

0.0458 

0.0963 


Table HI. Effect of angle of attack, sideslip, and elevons and flaperons on 

pitching-moment coefficient. 
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Table IV. Effect of angle of attack on stability and control parameters. 
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a, deg 

mg 

k = .135 

wmm 


mg 


20.8 

2.7173 

2.7262 

2.7169 

2.6812 

2.7182 


25.9 

2.0605 

2.3022 

2.4510 

2.4699 

2.6217 


30.8 

1.0483 

1.5909 

1.8777 

1.9985 

2.1274 

C L „ 

35.8 

0.4157 

1.0925 

1.2745 

1.4028 

1.6543 

a 

40.8 

0.3572 

0.7653 

0.9657 

1.0359 

1.2423 


45.9 

- 0.1074 

0.3471 

0.4584 

0.6013 

0.9278 


50.8 

- 0.2827 

0.0172 

0.2680 

0.2907 

0.4220 


55.9 

- 0.5163 

- 0.3175 

- 0.1771 

- 0.0864 

0.0863 


61.1 

- 0.7259 

- 0.6200 

- 0.4795 

- 0.4348 

- 0.2624 


20.8 

0.6274 

- 0.6749 

- 0.4411 

- 0.4391 

0.1617 


25.9 

6.8878 

2.8420 

1.3282 

0.5986 

0.4959 


30.8 

13.4790 

6.8415 

3.5802 

2.3319 

1.0400 

C Lq 

35.8 

17.0380 

7.8742 

4.5156 

3.2101 

1.4278 


40.8 

18.0300 

7.8265 

4.2304 

2.9332 

1.4082 


45.9 

10.8080 

5.3759 

3.5004 

2.2017 

1.2853 


50.8 

7.3799 

3.6589 

2.0937 

1.2329 

0.6319 


55.9 

4.0848 

2.6166 

1.5641 

1.0655 

0.2893 


61.1 

3.5036 

1.7856 

1.3260 

1.0968 

0.2086 


Table V. In-phase and out-of-phase components of lift coefficient. 
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20.8 

3.2712 

3.2778 

3.2707 

3.2352 

3.2971 


25.9 

2.8836 

3.1394 

3.3009 

3.3276 

3.5424 


30.8 

2.0690 

2.6814 

3.0051 

3.1475 

3.3558 


35.8 

1.5553 

2.3672 

2.5946 

2.7657 

3.1757 

Cjsj 

A ^a 

40.8 

1.5905 

2.1257 

2.3880 

2.5087 

2.8960 


45.9 

1.0970 

1.7409 

1.9091 

2.1241 

2.7033 


50.8 

1.0552 

1.5107 

1.9158 

1.9464 

2.3034 


55.9 

0.9679 

1.3127 

1.5457 

1.7166 

2.1429 


61.1 

0.9027 

1.1141 

1.4076 

1.5303 

2.0425 


20.8 

0.6466 

- 0.8186 

- 0.5106 

- 0.4904 

0.1744 


25.9 

7.4524 

2.9259 

1.3458 

0.5894 

0.5346 


30.8 

15.0792 

7.4842 

3.9208 

2.5651 

1.1490 

Cn 

35.8 

20.4228 

9.2855 

5.3797 

3.8608 

1.7763 


40.8 

23.1772 

9.9533 

5.4853 

3.8410 

1.9113 


45.9 

15.4526 

7.5281 

4.9940 

3.2200 

1.8646 


50.8 

11.8092 

5.7557 

3.3916 

2.1095 

1.2027 


55.9 

7.4302 

5.2961 

2.7569 

1.8374 

0.6297 


61.1 

8.0909 

4.2358 

2.3987 

2.0903 

0.3287 


Table VI. In-phase and out-of-phase components of normal-force coefficient. 
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a , deg 


mm 


in 


20.8 

0.5183 

0.5096 

0.5076 

0.5013 

0.4656 

25.9 

0.4529 

0.4775 

0.4937 

0.4941 

0.5023 

30.8 

0.1901 

0.2764 

0.3429 

0.3668 

0.3865 

35.8 

0.1090 

0.2378 

0.2325 

0.2419 

0.2278 

40.8 

0.1509 

0.1682 

0.1318 

0.1336 

0.0124 

45.9 

0.0323 

0.0706 

0.0253 

0.0458 

- 0.0816 

50.8 

- 0.0921 

- 0.0733 

0.0041 

0.0489 

0.0633 

55.9 

- 0.1755 

- 0.1386 

- 0.1546 

- 0.0612 

0.1958 

61.1 

- 0.2081 

- 0.1714 

- 0.1520 

- 0.1 195 

0.1636 

20.8 

- 0.8553 

- 1.1161 

- 1 .0749 

- 1 .0589 

- 1 .0527 

25.9 

- 0.1240 

- 0.7076 

- 0.8825 

- 0.9682 

- 1 .0761 

30.8 

1.7197 

0.4942 

- 0.1789 

- 0.261 1 

- 0.7195 

35.8 

2.3630 

0.5519 

- 0.3245 

- 0.4752 

- 0.7640 

40.8 

1.3537 

- 0.3968 

- 0.7955 

- 0.7878 

- 0.6640 

45.9 

- 0.1715 

- 0.5400 

- 0.4869 

- 0.6901 

- 0.6076 

50.8 

0.3532 

- 0.0418 

0.0527 

- 0.1078 

- 0.4099 

55.9 

0.0235 

0.1021 

0.1407 

0.1023 

- 0.2927 

61.1 

- 0.0310 

0.0538 

0.0567 

0.2266 

0.0463 


Table VII. In-phase and out-of-phase components of pitching-moment 

coefficient. 
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Measured 

Data 


Standard Error 


<*) 


mm 


max 


min 


5 ( c *. <°°>) 


max 


0.056 


0.074 


0.13 


0.19 


0.15 


0.21 




C N„ 


0.091 


0.12 


0.21 


0.30 


0.25 


0.34 


'm r 




0.034 


0.14 


0.069 

0.090 


0.17 

0.11 


* for index A a = L a ,L q ,N a ,N q ,m a ,or m q 


Table VIII. Minimum and maximum values of standard errors of estimated parameters. 


Measured 

Parameter 

Data 

*1 

bj,sec 1 

Tj ,sec 

Ct ? Ct 

17.2 ± 1.0* 

2.71 ±0.16 

0.368 + 0.023 

c Na ,c Nq 

17.1 ± 1.3 

2.73 + 0.21 

0.365 ± 0.028 

c c 

L m a ’ 

25.1 + 8.7 

1.86 + 0.64 

0.540 ± 0.190 


* the second value indicates the standard error 


Table IX. Estimated and computed parameters. 
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Parameter 

c 

'L 

C 

i 

'N 

e 

•w 

e 


e 0 

16.0 

1.4 


2.1 

01 

- 48.4 

4.8 

- 46.4 

7.1 

e 2 

33.6 

3.8 

32.4 

5.7 

0 3 

- 50.9 

11.3 

- 44.6 

16.9 

04 

- 2.0 

1.7 

- 2.6 

2.5 

05 

5.7 

3.9 

7.2 

5.8 

06 

- 3.4 

2.7 

- 4.1 

4.0 

07 

12.1 

0.51 

13.8 

0.75 

08 

- 46.2 

1.7 

- 51.8 

2.5 

0g 

35.7 

1.3 

38.9 

1.9 

010 

- 57.8 

3.4 

- 59.5 

5.0 


Table X. Estimated parmeters and their standard errors in Model II for measured lift and 

normal-force coefficients. 
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Figure 3. Effect of sideslip on longitudinal coefficients. 5 e =0 
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Figure 9. Concluded. 
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Figure 13. Amplitude and phase characteristics of C,/a transfer function 
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Figure 15. Variation of in-phase and out-of-phase components with sideslip angle 
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Figure 17. Estimated parameters of lift components. Model I 








Figure 18. Estimated parameters of normal-force components. Model I 
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Figure 21 . Measured and estimated in-phase and out-of-phase components of 
normal-force coefficient. Model I. 
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Figure 26. Estimated parameters of lift components 












Figure 28. Measured and estimated in-phase and out-of-phase components of lift 

coefficient. Model II. 
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Figure 29. Measured and estimated in-phase and out-of-phase components of normal-force 

coefficient. Model II. 
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Figure 30. Measured and predicted in-phase and out-of-phase components of lift 

coefficient. Model II, k=.190. 
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